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The synthesis of ultra-long lived acoustic phonons in a variety of materials and device geometries 
could enable a range of new coherent information processing and sensing technologies; many forms 
of phonon dissipation pose a barrier to this goal. We explore linear and nonlinear contributions 
to phonon dissipation in silica at cryogenic temperatures using fiber-optic structures that tightly 
confine both photons and phonons to the fiber-optic core. When immersed in helium, this fiber 
system supports nearly perfect guidance of 9 GHz acoustic phonons; strong electrostrictively me¬ 
diated photon-phonon coupling (or guided-wave stimulated Brillouin scattering) permits a flexible 
form of laser-based phonon spectroscopy. Through linear and nonlinear phonon spectroscopy, we 
isolate the effects of disorder-induced two-level tunneling states as a source of phononic dissipation 
in this system. We show that an ensemble of such two-level tunneling states can be driven into 
transparency-virtually eliminating this source of phonon dissipation over a broad range of frequen¬ 
cies. Experimental studies of phononic self-frequency saturation show excellent agreement with a 
theoretical model accounting for the phonon coupling to an ensemble of two-level tunneling states. 
Extending these results, we demonstrate a general approach to suppress dissipation produced by 
two-level tunneling states via cross-saturation, where the lifetime of a phonons at one frequency 
can be extended by the presence of a high intensity acoustic beam at another frequency. Our mod¬ 
eling and measurements suggest that Rayleigh scattering dominates phonon losses for the longest 
lifetimes achieved in our system. Although these studies were carried out in silica, our findings are 
quite general, and can be applied to a range of materials systems and device geometries. 


I. INTRODUCTION 

Access to new regimes of classical and quantum dy¬ 
namics hinge upon our ability to create and manipu¬ 
late ultra-long lived coherent excitations in electromag¬ 
netic, optical, and phononic domains DMG2- In partic¬ 
ular, ultra long-lived phonon modes have been identi¬ 
fied as a crucial new resource by quantum information, 
optomechanics, and precision metrology communities 
000IHHI9]. To this end, a variety of systems, ranging 
from nano- and micro-scale phononic devices to resonator 
technologies of centimeter-scale, have harnessed remark¬ 
able phonon coherence times [21131131 IT21 121)1123] : however, 
radical improvements in performance are possible if tech¬ 
nical and fundamental sources of dissipation are mastered 
0 0 EH TB- This realization has spawned a resur¬ 
gence of interest in the fundamental origins of phonon 
dissipation at cryogenic temperatures [22] [24H29] . 

A ubiquitous source of dissipation arises from disorder- 
induced defects. Some of such defects possess quantized 
energy spectra, and can exchange energy with electro¬ 
magnetic, optical, and phononic fields [22J :23] I30U35] . 
When phonon-active, such defect states can absorb and 
emit phonons just as atoms absorb and emit light. Dis¬ 
sipation by such defects pose a fundamental limit to 
phonon lifetimes in acoustic media at cryogenic temper¬ 
atures. Phonon-active defect states have been exten¬ 
sively studied in (highly disordered) amorphous media 
[20113 03 [MSa; however, their deleterious effects also 
appear within highly ordered crystalline systems and in 
systems with material interfaces 0 eh [ 22 ] . This form 
of dissipation bars access to new regimes of classical and 
quantum dynamics, central to a range of emerging tech¬ 


nologies 0 ;TU]. 

In this work, we examine the dynamics of an ensem¬ 
ble of phonon-active defects using a guided wave geom¬ 
etry that produces tight confinement of both light and 
acoustic waves (Fig 1). Strong photon-phonon coupling 
within this system permits noninvasive optical excita¬ 
tion and interrogation of high frequency (~ 9.2 GHz) 
phonons. Tight confinement of acoustic modes pro¬ 
duce high phonon intensities (600 W/m 2 ) in the core of 
this waveguide, permitting frequency selective nonlinear 
phonon spectroscopy with modest (~mW) optical pow¬ 
ers over a range of temperatures (1.1-300 K). Building on 
established models, we elucidate the nature of phonon- 
defect interactions in our system, allowing us to extract 
defect density, coupling strength, and a range of other pa¬ 
rameters that capture the dynamics of the defect ensem¬ 
ble. We show that this ensemble of defects can be driven 
into transparency in the strong-field limit yielding an es¬ 
timated factor 45 suppression of defect-induced dissipa¬ 
tion at 1.1 Kelvin. In this limit, defect-induced dissipa¬ 
tion is a negligible source of loss within our guided-wave 
system; remarkably high phononic Q-factors (>12,000) 
and decay lengths (>lmm) are achieved. Building on 
these findings, material engineering and more sophisti¬ 
cated schemes could provide a path toward radically en¬ 
hanced coherence times [33] in a range of systems, as the 
basis for emerging quantum information technologies. 

Phonon dissipation can be suppressed by lowering the 
system temperature (below 150K in our system), and 
may in itself be a viable strategy to achieve low phonon 
losses in crystalline material which have low defect den¬ 
sities. However, for amorphous (and even crystalline) 
materials such as silica or silicon nitride, a tempera- 
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ture will be reached below which the acoustic damp¬ 
ing will cease to decrease as the temperature is lowered 
further ISHEl S5H13- This effect is believed to arise 
from resonant absorption by hypothetical two-level tun¬ 
neling states (TLSs) [35HH]; the dominant contribution 
to phonon dissipation in many materials at low temper¬ 
atures. 

Tunneling states are hypothesized to arise from a sub¬ 
set of atoms that inhabit asymmetric double-well poten¬ 
tials (see Fig[l]). A perturbation of this potential 

by strain enables phonon absorption, and at low tem¬ 
peratures when these TLSs condense into their ground 
states (see Figs. [I] and [2]) they are capable of resonant 
absorption of phonons with energy E. The two-level na¬ 
ture of this decay channel leads to non-linearity in the 
phonon dynamics [3BJ f55H45j and provides a way to break 
through the dissipation floor established by resonant ab¬ 
sorption. 

We demonstrate that by working with high phonon 
intensities resonant absorption can be saturated, signif¬ 
icantly extending phonon lifetimes. This self-frequency 
saturation is accomplished by increasing the amplitude 
of an acoustic mode (at a single frequency) to drive the 
TLSs into transparency. 

Two-level tunneling states have generated much con¬ 
temporary interest in their own right for the possibil¬ 
ity of beneficial implications. In particular, it has been 
shown that they can engender the dynamics of nano¬ 
electromechanical systems with nonlinearity, even in the 
single phonon regime |44j ■ and their relatively long coher¬ 
ence times have facilitated their use as a quantum mem¬ 
ory 134) . In light of this recent interest, our exploration 
of TLSs and their manipulation may provide valuable 
information for a variety fields. 


II. OVERVIEW 

A. Stimulated Brillouin scattering 

Our system is a 2.2 cm segment of ge-doped Nufern 
UHNA-3 optical fiber. This fiber’s high germanium con¬ 
centration (~ 44 wt%) guides light exceptionally and 
produces a large longitudinal sound speed contrast be¬ 
tween core UL,core = 4,740 ± 68 m/s [48] and cladding 
*’L,ciad = 5,944 m/s [45] that promotes acoustic guid¬ 
ance. The strong confinement and guidance of optical 
and acoustic modes (see Figs. [l]&[3| allows for the effi¬ 
cient excitation of phonons via stimulated Brillouin scat¬ 
tering (SBS) and which provides a means to characterize 
our system [55M57] , 

Stimulated Brillouin scattering is a resonant process 
involving the interaction of light and sound. It occurs 
in media exhibiting photoelasticity and its complement 
electrostriction which are respectively characterized by a 
change in refractive index upon strain and by an optical 
intensity-induced elastic deformation. Hence, the beat 
tone formed by two counter propagating optical fields can 


drive acoustic waves through electrostriction, and once 
created, acoustic waves act as moving Bragg gratings 
which reflect and Doppler shift incident optical beams 
as a consequence of photoelasticity. 

The concert of these two effects transfer energy be¬ 
tween the two participating pump and Stokes optical 
fields, and the acoustic field with respective angular fre¬ 
quencies uip, to s, Cl and wavevector magnitudes k p , kg, q. 
The transfer of energy is only significant when the three 
fields are phase matched (see Fig. [l]g-h); for backward 
SBS in ID ) 5Sj 

UJp = ws + fi (1) 

k p = q-k s , (2) 

alternatively the phase-matching conditions can be 
viewed as energy and ‘momentum’ conservation. The 
minus sign in front of kg in the second equation above 
indicates that the Stokes photon and the pump counter- 
propagate (the phonon and the pump copropagate). For 
a given pump frequency we can approximate the phonon 
frequency, excited via SBS, by assuming linear dispersion 
for light (to = ( c/n)k ) and sound (Cl = vq) 


Cl « 


2 nv 

-w p 
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(3) 


where n is the effective index of refraction, v is the modal 
sound speed, and c is the speed of light. 

In our system SBS results in the transfer of energy 
from the pump to Stokes beam. Net Stokes amplifica¬ 
tion can be measured by balanced detection. For such 
measurements it is prudent to work in the weak signal 
regime, given by g B P p L <C 1, where g b is the Brillouin 
gain (gs ~ O.G(Wm)^ 1 at room temperature in UHNA-3 
fiber), P p is the pump power, and L is the length of the 
fiber under test (FUT). In this case the power transferred 
to the Stokes beam APg ? which is directly acquired by 
balanced detection, is small and is given by 


(T/2) 2 g B P p P s L 
(Cl-tui M ) 2 + (r/2) 2 


(4) 


where we have dropped a negligible correction arising 
from the optical loss, Ps is the input power of the Stokes 
beam, the angular frequency wim is the angular frequency 
detuning between the pump and Stokes beams which is 
experimentally set by an intensity modulator (IM) see 
Fig. 0 and 1/T is the phonon lifetime. By sweeping 
wim balanced detection supplies the peak Brillouin gain, 
and the phonon frequency and lifetime which characterize 
SBS’s Lorentzian frequency response. 


B. Phonon dissipation in glasses 

Dissipation of phonons in glasses has several origins: 
from multi-phonon interactions to scattering by defects. 
Among the latter of these are hypothesized TLSs that at¬ 
tenuate acoustic waves through the processes of resonant 
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FIG. 1. a) Illustration of the microscopic structure of silica glass and some candidate defects proposed as the origin of two- 
level tunneling states (l)-(3). Individual TLSs are characterized by an asymmetric double-well potential for the generalized 
coordinate x with well-separation d, asymmetry A, barrier height V, and the average oscillation frequency of the two individual 
wells uj c . Tunneling between the two-wells is characterized by the tunneling strength Ao (defined below). Excited |e) and 
ground \g) energy eigenstates of an uncoupled TLS are split by energy E = \J A 2 + A ( 2 ,. b) Guided optical and acoustic beams 
used for the stimulated Brillouin spectroscopy of our system, c) & e) define the material properties and geometry of the fiber 
system under study. LHc stands for liquid helium. Optical d) and acoustic f) intensity profile for the modes participating in 
Brillouin scattering. MFD is mode field diameter, g) spatial and h) temporal phase matching conditions between the optical 
and acoustic waves in SBS (see text). 


scattering and relaxation absorption. To understand the 
origins of these two mechanisms we lay out the tunneling 
state model and discuss its consequences for the dynam¬ 
ics of phonons. For those seeking more details, complete 
derivations of all results in this section can be found in 
the Appendix. 

A tunneling state is characterized by a double-well po¬ 
tential of asymmetry A and the overlap energy Ao = 
hu> c e~ x where Hu> c is roughly the average zero point en¬ 
ergy of the two wells and A = '/ 2 ™ vd characterizes the 
extent of wave function overlap between the two wells 
with to the mass of the atom(s) comprising the TLS, V 
the barrier height, d the ‘distance’ between the double- 
well’s minima (see Fig. IT]) , and where h is Planck’s con¬ 
stant divided by 2ir. Given a finite barrier height the 
atom(s) may tunnel between the two minima of the po¬ 
tential. The parameters A and A are assumed to be uni¬ 
formly distributed over the ensemble of TLSs resulting in 
a constant density of states G(A, A) = P, or in terms of 
the overlap energy /(A,Ao) = P/ Ao- This assumption 
of uniformity correctly predicts the linear in temperature 
behavior of the specific heat and the anomalous thermal 
conductivity of glasses at low temperatures loilffTA] . 

At low temperatures the TLS+phonon system, includ¬ 
ing the strain-induced perturbation of the asymmetry, 
can be modeled with the following Hamiltonian 


H = H, 


P h 


E 
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2 (Ej + Dj ' £( r j)) a z,j + Mj ■ £(r j)cr x j 


( 5 ) 


[42] where the physics of the jth TLS is approximated 
as an effective two-level system using the Pauli matri¬ 
ces cTij with i = {x,y, z}. The elastic strain field of 
polarization q is labeled as £ r /, the sum over j counts all 
TLSs in the glass at various positions r j with energy split¬ 
ting Ej and coupling parameters D v j = 2(Aj/Ej)'y r) and 
M v j = (Ao j/Ej)^, Pph is the free Hamiltonian for the 
phonons (see Appendix), and the ‘dot product’ Dj ■£ ab¬ 
stractly represents a sum over polarizations ^ D Vi j£ v . 
Generally, the deformation potential tensor, quantifying 
the TLS-strain coupling, is unique to each TLS and ori¬ 
entation dependent; we have ignored this complication 
above where the TLS-phonon interaction is character¬ 
ized by a polarization dependent deformation potential 
constant 7 ,,. 


C. Resonant phonon absorption by TLSs 


Tunneling states can resonantly interact with phonons 
though three processes: stimulated absorption, and spon¬ 
taneous and stimulated emission. The relative magnitude 
of each of these processes determines the phonon dissipa¬ 
tion rate and depends on temperature and the intensity 
of the acoustic field. 
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a) High temperature transparency 


Pe/Pg « 1 



b) Low temperature absorption 


Pe/Pg ~ 0 




Radius (pm) 

FIG. 3. Normalized fundamental mode acoustic and optical 
intensity, and index profile in UHNA3 fiber. Solid lines are 
obtained from finite element simulations, and red dots from 
analytical calculations of the acoustic modes. Inset: Cross 
sectional view of the acoustic and optical intensity. The small 
core size and good opto-acoustic mode overlap ensures high 
optical and acoustic intensities as well strong optomechanical 
coupling. 


c) Saturation-induced transparency 
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FIG. 2. Resonant absorption of phonons of angular frequency 
Q by TLSs (illustrated by green spheres), a) At high tem¬ 
peratures (ksT 3> lift) the TLSs addressed by the acoustic 
wave have almost equal probabilities of being in the excited 
or ground states, and hence resonant absorption is compen¬ 
sated by stimulated phonon emission, b) At low tempera¬ 
tures (ksT <C hQ) the TLSs addressed by the acoustic wave 
condense to their ground states. At sufficiently low acoustic 
intensities TLSs, having resonantly absorbed a phonon, de¬ 
cay via spontaneous emission thus attenuating the incident 
acoustic wave in the process, c) At low temperatures and 
high phonon intensities a given TLS is interrogated by mul¬ 
tiple phonons during its excited state lifetime. Hence, an ex¬ 
cited TLS may decay via stimulated phonon emission which 
preserves the coherence of the incident phonon. At high in¬ 
tensities resonant absorption is compensated by stimulated 
emission. 


1. Weak fields 

First, we consider weak fields which is roughly charac¬ 
terized by a mean free time between TLS-phonon inter¬ 
actions that is long compared to the excited state lifetime 
of the TLS, and which will be defined quantitatively be¬ 
low. For such low intensity acoustic waves of angular 
frequency Q the golden rule can be used to find the TLS- 
induced decay rate (see Appendix) 




pv 
TlP^Q 
pv 2 


( 6 ) 


(i) («) 


where the polarization dependence of the defor¬ 
mation potential 7 and the sound speed v have been 
suppressed, ks is the Boltzmann constant, and p is the 
material density. F res is characterized by two regimes 
occurring at high and low temperatures T. This is elu¬ 
cidated by the temperature dependence of the phonon 
decay rate on the thermal equilibrium population inver¬ 
sion P e —P g of the TLSs at energy Ml and temperature T, 
where P e and P g are the probabilities to find a TLS in the 
excited and ground state, respectively. At low tempera¬ 
tures Ml ksT the decay rate is maximized since the 
TLSs are found entirely in their ground state (P g —> 1), 
and thus (i) stimulated absorption dominates. In the low 
intensity regime TLSs that have absorbed phonons decay 
through spontaneous emission which attenuates coherent 
phonon beams. As the temperature is raised the prob- 
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ability to find the TLS in its excited state grows which 
opens the possibility for (ii) stimulated phonon emission 
which coherently amplifies the sound amplitude. At high 
temperatures these two processes compensate each other 
exactly and resonant absorption is suppressed. 


2. Strong fields 

At high acoustic intensities perturbation theory is no 
longer adequate to calculate the phonon decay rate. In 
this regime one must solve the coupled Heisenberg equa¬ 
tions of motion for the TLSs and the acoustic field (see 
the Appendix for a complete derivation). The backreac- 
tion of the TLSs on the phonons is characterized by a 
complex susceptibility that modifies the phonon dynam¬ 
ics in two ways: by resonant absorption induced dissipa¬ 
tion and a frequency shift. At high acoustic intensities 
the dissipation rate due to resonant absorption is given 
by 


F res _ nP"/ 2 n tanh(hfl/2k B T) ^ 

pv 2 y/1 + J/Jc ' 1 

[42] where J is the acoustic intensity, and J c is the critical 
intensity that demarcates the boundary between weak 
and strong fields. The critical intensity 


h 2 pv 3 
2 7 2 TiT 2 


( 8 ) 


is characterized by the time scales Tj and T 2 which 
are phenomenological upper state lifetime and dephasing 
times for the TLSs with energy E ss hQ, (see Appendix). 

As with the case of weak fields, this result for the 
phonon dissipation can be interpreted as a competition 
between stimulated emission and absorption. In the 
strong field case however the effective population inver¬ 
sion for the TLSs which interact with a phonon of fre¬ 
quency fl is given by the nonequilibrium steady-state 
value 


tanh(M2/2fc#T) 
\/l + J/ Jc 


(9) 


At high intensities J/J c > 1 a TLS, excited through res¬ 
onant absorption, will encounter several phonons within 
its upper state lifetime. Thus, stimulated emission can 
become a dominant decay channel and resonant absorp¬ 
tion will be suppressed. 

Resonant absorption also leads to a (acoustic intensity- 
insensitive) frequency shift given by 


Aw res (T)—Aw res (r 0 ) 


P 7 2 H 

pv 2 


In 


- ReT 


/1 nn \ 

\2 + 2t rik B T J 


/ m 

[k^T 

-(T^T 0 ), (10) 


where T is the digamma function i22j . 


D. Relaxation absorption 


In addition to resonant absorption TLSs can atten¬ 
uate coherent acoustic waves via relaxation absorption. 
This process results from a modulation of the TLS’s en¬ 
ergy splitting by the presence of a time-dependent strain 
field, i.e. Ej(t) —> Ej + D v ^£ iV {t). As Ej(t) changes, the 
TLS can equilibrate with the surroundings by absorbing 
or releasing energy. Thus, the instantaneous ‘equilib¬ 
rium’ population inversion is modulated in time. For our 
system (see Appendix) relaxation-absorption results in 
dissipation of phonons characterized by the decay rate 


prel 


7T 3 Pit 
24 p 2 v 2 h 4 



C k B T ) 3 , 


( 11 ) 


and a temperature-independent frequency shift. 


III. EXPERIMENTAL SETUP 

The system is characterized by pump-probe measure¬ 
ments performed using balanced detection and lock-in 
amplification. Our apparatus consists of a 1548.963 nrn 
source that is split into two optical lines (see Fig. [4|; 
one reserved to act as a pump, and the second to act as 
a probe. The pump line is subsequently amplified and 
pump and probe polarizations are aligned. The pump 
beam is then power modulated at the fixed frequency 
Umod f° r lock-in detection before being sent through the 
FUT. The probe beam is sent through an intensity mod¬ 
ulator which upon exit is filtered to isolate a single side¬ 
band. To achieve common mode noise rejection using bal¬ 
anced detection, the probe beam is split into two arms. 
One arm passes probe light through the FUT and is am¬ 
plified via SBS whereas the other arm acts as a reference. 
A variable attenuator on the reference arm is adjusted so 
that both probe arms have the same optical power in the 
absence of gain via SBS. The difference in power of the 
balanced probe arms yields the net gain experienced by 
the probe beam passing through the sample. 

Modulation of the pump power generates a fixed- 
frequency side-band of the amplified probe which is de¬ 
tected with a lock-in amplifier. To measure Brillouin gain 
spectra (BGS) the modulation frequency wim is swept 
across the Brillouin resonance generating a Lorentzian 
response. 

The temperature of the fiber is controlled using a 
double walled cryostat, and a large copper heat sink. 
The samples are mounted in a shallow mail slot passing 
through the block. The block size was chosen to ensure 
temperature stability, and to allow a minimum of 10 Bril¬ 
louin lineshape measurements per 100 mK as the block 
warms slowly to room temperature. The system is cooled 
using liquid Helium and evaporative cooling resulting in 
lowest achievable temperatures in the neighborhood of 1 
K. 
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FIG. 4. Experimental setup to measure SBS with balanced 
detection. The acronyms above are defined as: EDFA Erbium 
doped fiber amplifier, VOA variable optical attenuator, RF 
radio frequency signal generator, BD balance detector, and 
FBG fiber Bragg grating. We use a DFB (distributed feed¬ 
back) laser with a wavelength of 1548.963 nm and a linewidth 
of 5 kHz. Our probe line is synthesized by sending pump 
laser light through an IM modulated in the neighborhood of 
9.2 GHz, the carrier and the high frequency side-band are 
removed by operating at the null point bias point of the IM 
and band-pass filtering. Probe light is split in two arms, one 
acting as a reference and the other traversing the FUT. The 
power modulated pump results in a time harmonic amplifi¬ 
cation of the Stokes beam at f! mo d which is acquired with 
lock-in detection. 


To study the saturation of losses arising from the TLSs 
the intensity of the sound field inside the FUT was swept 
at fixed temperature. In the weak signal regime the 
steady-state phonon intensity follows the powers of the 
optical fields 


J(wim) ~ w —— w—G b(wim)-P p -Ps, (12) 

r uj s A eff 


where A B s ~ 1 . 6 /im 2 (estimated from simulations, see 
Fig. § is the phonon mode area. Hence, the phonon 
intensity can be controlled by changing the power of the 
optical driving fields. By sequential adjustment of the 
variable optical attenuators on pump and probe lines the 
phonon intensity generated through SBS can be swept 
through more than 4 decades of dynamics range. The 
maximum pump and probe powers were used to deter¬ 
mine the 2.2 cm length of the FUT; this length en¬ 
sures that all measurements were performed in the weak- 
signal regime. This fiber length allowed access to high 
phonon intensities while simultaneously preventing in¬ 
homogeneous broadening from non-linearity induced by 
strong backward scattering via SBS. 


IV. RESULTS 

A. Temperature and Intensity Dependence of 
Phonon Losses 

Brillouin gain spectra corresponding to the fundamen¬ 
tal acoustic mode in a 2.2 cm segment of UHNA-3 fiber 
were continuously acquired as the system slowly warmed 
to room temperature from 4.17 K. Pump and probe pow¬ 
ers were set to 35 mW and 0.550 mW, respectively to en¬ 
sure SBS measurements in the weak signal regime. The 
BGS were binned in 100 mK steps, averaged, and fit to 
the Lorentzian model given by Eq. [4j This analysis pro¬ 
vides the phonon frequency fI and the dissipation rate 
r as a function of temperature. The red data points in 
Figure [5|i show the dissipation rate of the fundamental 
acoustic mode as a function of temperature. Three rep¬ 
resentative BGS are shown in Fig. B- To judge the 
relative importance of resonant absorption by TLSs we 
have plotted the ratio P e /P g as a red curve (top). 

For low phonon intensities (J <C J c ) the linewidth be¬ 
gins to level off below 4K (red points in Figs, 
and then begins to increase as the temperature is lowered. 
This gray region in Figs. B &j5jr, where P e /P g < 0.9, in¬ 
dicates the temperature range where resonant absorption 
by TLSs begins to dominate the acoustic damping. 

After cooling the system to 1.1 K, we let the fiber sys¬ 
tem slowly warm up to 4.17 K. For each 100 mK rise 
in temperature we acquire BGS as a function of phonon 
intensity. The data is analyzed by binning and averaging 
as described above, and results in dissipation rate as a 
function of both temperature and phonon intensity. This 
data is presented in [5 Jd & B-B shows phonon dissipa¬ 
tion rate as a function of temperature. The three sets of 
data correspond with three distinct settings of the opti¬ 
cal powers that generate phonons of low, moderate, and 
high intensities (as compared to J c ). Normalized BGS 
for the three settings is shown in Fig. [5}i for the lowest 
temperature of 1.1 K. Figure[5j: shows phonon dissipation 
rate as a function of phonon intensity for three different 
temperatures (black, orange and tan points). It is qual¬ 
itatively clear from Figs. that at temperatures 

lower than 4 K phonon dissipation is suppressed at high 
phonon intensity; this indicates saturation of the reso¬ 
nant absorption due to TLSs and is made evident by the 
theoretical plot of P e / P g (at top), computed using Eq. 

El 

In Figs. B f, & g we compare our measurements with 
the tunneling state model. In Fig. B a model of the 
phonon dissipation rate given by 

r = r res + r 0 (13) 

is fit to the data where T 0 is an offset parameter that 
represents all intensity independent background acoustic 
losses including relaxation absorption, Rayleigh scatter¬ 
ing, phonon-phonon interactions, etc. Three parameters 
are employed to fit the data: 1 ) P 7 2 ( 7 l being the de¬ 
formation potential for longitudinal waves), 2) J c , and 
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FIG. 5. a) Dissipation rate of the fundamental acoustic mode at low intensity as a function of temperature, b) Phonon 
dissipation rate as a function of temperatures below 6 K. The three data sets below 4 K correspond with different pump 
and probe power settings corresponding with low, moderate, and high phonon intensities, as compared to J c . c) Phonon 
dissipation rate as a function of acoustic intensity. Each curve has a fixed temperature, d) Brillouin gain spectra as a function 
of temperature, e) Normalized Brillouin gain spectra as a function of intensity at 1.1 K. f) Comparison of the measured 
Brillouin frequency vs. temperature with the theoretical prediction given by Eq. 10 g) Critical intensity J c and background 


losses Fq as a function of temperature. The parameters for the theory can be found in Tab.l 


3) the offset r 0 . The parameters P and 7 ^ are assumed 
to be constant. In distinction, J c and To are expected 
to depend upon temperature. In particular, J c scales in¬ 
versely with the effective decay rates TjT 2 . The upper 
state lifetime Tj can be approximated with the Fermi’s 
golden rule (see Eq. Cl of the Appendix), but an esti¬ 
mation of the temperature dependence of T 2 is beyond 
the scope of this work. Instead we fit J c (gray dots in 
Fig [HJr) to a power law J c = aT b that is shown as a solid 
gray line in In Fig. [5f the measured shift in Bril¬ 
louin frequency (gray points) referenced to To = 1.09 K 
is compared to Eq. [To] The fitted values obtained from 
Fig. [5J: were used as inputs. The offset r 0 is plotted as 


blue points in Figs, [h]* and contains all background losses 
which we model as 

r 0 = r rel + r BG . (14) 

where T BG is a constant representing the background dis¬ 
sipation floor for our system. The model is fit to the data 
by adjusting the offset and the value of the longitudinal 
deformation potential 7 ^ (we assume that the transverse 
deformation potential is given by 7 ^ ~ 7^/2 [35]). P"l\ 
is held fixed to the value obtained from the data analysis 
in Fig. [5jc. 

The model parameter values obtained from this analy¬ 
sis are tabulated in Tab. [l]and their values are compared 



























































44% wt. ge-doped silica 

vitreous silica 

p 

2, 666 kg m -3 

2, 202 kg m" 3 US] 

VL 

4, 760 ± 68 m s" 1 [SI 

5,944 m s 1 [49| 

a 

0.9 W m" 2 KT ft 


b 

2.6 


Jc 

1.2 W m~' 2 


Pit 

1.6 x 10 Y J in -3 

1.3 x W J m" 3 [5S] 

p 

23 

6.85 [58] 

7 L 

0.5 eV 

0.86 eV [58] 

VTiT 2 

10 ns (est. with Eq 81 

*53 //s [37] 

T 

« 79 ns (est. with Eq Cl 1 

*200 ns (37j 

t 2 

« 1.3 ns 

*14 /xs [37] 


TABLE I. Comparison of the properties of ge-doped silica and 
silica. All temperature and frequency dependent quantities 
are listed for T = 1.1 K and Q = (2-7r)9.188 GHz. The product 
T 1 T 2 is extracted from J c defined by Eq. [8] The parameter 
T\ above is calculated from Eq. |C1| and then subsequently 
employed to extract T 2 from J c . P is in units of 10 44 J -1 m -3 . 
*The values of Ti and T 2 for vitreous silica were measured at 
T = 20 mK and for frequency 0.68 GHz. Hence, these values 
should not be directly compared the values listed for ge-doped 
silica. Using Eq. |Cl| and the fitted tunneling state parameters 
Ti at 20 mK and for 0.68 GHz phonons is estimated at 665 
/rs in UHNA-3, this compares well with the listed value for 
vitreous silica. 


with those of vitreous silica. There are several remarks 
in order. The data of Fig. is well-described by the 
tunneling state model, and the computed Brillouin fre¬ 
quency shift using the fitted value of from Fig. [Hj; 
compares well with our measurements, however the the¬ 
ory and experiment begin to diverge at higher tempera¬ 
tures. Relaxation absorption accounts well for the inten¬ 
sity independent background, and the fitted value for 
is comparable to silica. The fitted parameters and their 
theoretical relationships can also be used to estimate the 
density of TLSs P and the dephasing rate T 2 . 

As a final note, the values of Ti and T 2 should not be 
directly compared with those of vitreous silica listed in 
Tab. [T] The phonon frequencies and the temperatures 
were distinct from our measurements; for an appropriate 
comparison Ti can be estimated from the minimum TLS 
upp er st ate lifetime computed using perturbation theory 
Eq. Cl with Aq —> E, taking T = 20 mK and tt/2n to 
be 0.68 GHz, and using the fitted TLS parameters for ge- 
doped silica. We find Ti a; 665 /is which is comparable 
to the reported values for vitreous silica in Tab. |T] 


V. DISCUSSION 

Acoustic dissipation imposed by defects presents a bar¬ 
rier to attaining low loss acoustic modes that are criti¬ 
cal to a range of technologies. To overcome/understand 
these challenges we explored defect-induced dissipation 
in a guided wave system using nonlinear phonon spec¬ 
troscopy. Theory-experiment comparison using the tun¬ 


neling state model quantifies the influence of defects in 
our system. We found that the large phonon intensities, 
made possible by the tight acoustic confinement, permit 
access to the TLS-induced regime of nonlinear phonon 
dynamics. We demonstrated self-frequency saturation, 
where TLS resonant absorption is driven to transparency 
as the intensity of the acoustic beam becomes large. 

At the highest intensities we encounter a dissipation 
floor To due to relaxation absorption T rel and a roughly 
temperature independent offset Tbg- At the lowest tem¬ 
peratures and highest intensities we estimate that relax¬ 
ation absorption makes up only 6% of the dissipation rate 
floor, and thus we believe we have reached the neighbor¬ 
hood of the smallest phonon dissipation rates possible 
in UHNA3 fiber estimated to be (27 t) 650 kHz.There are 
several possible explanations for Teg such as; phonon- 
phonon scattering, anchor losses, inhomogeneous broad¬ 
ening due to irregularities in the fiber geometry along its 
length, acoustic leakage due to imperfect guiding, and 
Rayleigh scattering due to disorder and scattering by 
dopants. Phonon-phonon scattering can be estimated 
from the Landau-Rumer theory using Fig. 2 of [9] and 
gives negligible dissipation, and anchor losses are likely 
unimportant as the mode is highly confined to the core. 
Large amounts of inhomogeneous broadening are unlikely 
given the relatively short length of the fiber and the ob¬ 
served Lorentzian shape of the BGS, and acoustic leakage 
was estimated in simulations to be too small. 

We attribute the dissipation floor to Rayleigh scatter¬ 
ing. There are several pieces of corroborating evidence 
for this hypothesis; the fiber is doped with a large con¬ 
centration of germanium that will produce large stochas¬ 
tic variations of density and sound speed in the fiber 
core, Fbg is insensitive to temperature and intensity, and 
an estimation of Rayleigh scattering in Top High Qual¬ 
ity alpha synthetic quartz for 9.2 GHz phonons gives 
a lower bound for the Rayleigh scattering dissipation 
rate (27 t) 200 kHz. By working at lower frequencies (as 
Rayleigh scattering scales with the fourth power of the 
frequency) or with fibers with lower defect densities much 
longer phonon lifetimes may be possible. 

Acknowledgements: Primary support for this work 
provided by NSF grant DMR 1119826. 


Appendix A: Acoustic guidance in fiber 


The high germanium concentration (~ 44 wt%) in 
UHNA-3 fiber produces a large longitudinal sound speed 
contrast between core fL, C ore = 4, 740 ± 68 m/s |3S] and 
cladding unclad = 5,944 m/s jiS] that promotes acoustic 
guidance. However, ideal guidance is not expected since 
the shear wave velocity in the cladding Us jC iad = 3, 764 
m/s [J5j is smaller than UL,core> and hence leads to mode 
conversion at the core-clad interface [50] that leads to 
acoustic leakage. 

To understand the limits of acoustic guidance we stud¬ 
ied the axial-radial acoustic modes of our system with full 
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vectorial simulations and semi-analytical calculations. 
We solved for acoustic eigenfrequencies of a core region 
(diameter 1.8 /im) embedded in a finite cladding (diam¬ 
eter 125 (im) surrounded by liquid helium. A perfectly 
absorbing boundary condition was implemented adiabat- 
ically in the liquid region at large separation from the 
core to prevent spurious reflections from the simulation 
boundaries. We employed a shear velocity in the core 
region us,core ~ 3092 m/s which was estimated by inter¬ 
polating between the values of pure silica and germania, 
and the liquid region was ascribed the acoustic proper¬ 
ties of liquid helium. The core, cladding and liquid were 
assumed lossless other than the absorbing layer used to 
model mode leakage. 

Acoustic energy that enters the liquid region leaves 
the system irreversibly and hence reduces the acoustic 
guidance. The signature of this leakage is an imaginary 
part to the eigenfrequencies that quantifies the dissipa¬ 
tion rate. Our simulations correctly predict the acoustic 
frequency and yield negligible leakage (much smaller than 
the dissipation rate observed in experiment). Hence, we 
conclude that the dissipation of phonons in our system is 
due almost entirely to internal friction of the fiber as op¬ 
posed to leakage. In this sense we describe our phonons 
as well-guided. 

It should be emphasized that acoustic guiding in fiber 
is not guaranteed. Germanium doping has the fortunate 
consequence that increased doping leads to higher index 
of refraction while simultaneously decreasing sound speed 
m EH, thus enhancing the confinement of light and 
sound simultaneously. Titanium dioxide and phospho¬ 
rous pentoxide are other dopants which simultaneously 
enhance guiding of light and sound, while boron triox¬ 
ide and fluorine doping enhance the confinement of light 
while raising the sound speed [52]. Fibers doped with the 
former will possess leaky acoustic modes. At room tem¬ 
perature, where phonons are strongly damped in silica, 
one may posit that the guiding of sound in fiber plays 
a minor role as the mean free path is on the order of 
the acoustic mode field diameter. However, the existence 
of higher order guided phonon modes can be observed 
at room temperature [531 154] which shows that acous¬ 
tic modes in fiber cannot be completely understood in 
terms of bulk properties without considering boundaries. 
We have also simulated axial-radial modes in a system 
composed of a core embedded in an infinite cladding re¬ 
gion. In this case the computed acoustic leakage rate 
exceeds the observed phonon dissipation rate suggesting 
that the cladding liquid interface plays an important role 
in the observed acoustic guidance. In light of this fact it 
is important note that the leakage rate we derived from 
simulation for the finite fiber case is only a lower bound as 
we have not accounted for nonuniformity in the geometry 
of the fiber over its length. Irregularities in cladding and 
core diameters will contribute to inhomogenous broad¬ 
ening that will manifest as an effective increase in the 
leakage rate. 


Appendix B: Tunneling state theory of phonon 
dissipation in glasses 


In this section we summarize the underpinnings of the 
tunneling state theory. 

At low temperatures fc^T <C Hco c a TLS can be effec¬ 
tively represented by the Hamiltonian 


Htls 


1/A A 0 \ 
2 \A 0 -A) 


(Bl) 


where the states (L\ = (0,1) and (R\ = (1,0) correspond 
with position states with the particle localized in the left 
(‘L’) or right (‘f?’) well and which we approximate as 
being orthogonal. Diagonalizing the Hamiltonian above 
gives the corresponding energy eigenstates 


°) = 


An 


y/2E(E + A) . 

Ap 

V 2E ( E - A ) 


|L) + | R) 

Ao 

I L) + | R) 


(B2) 


which reveal that the stationary states of a free TLS are 
spatial superpositions between the wells and have ener¬ 
gies ±\E = ±^\/A 2 + Ag, see Fig. |TJ >. Such a TLS is 
perturbed by a local strain; the dominant effect of this is 
a shift in the asymmetry A that p ertu rbs the diagonal el- 


B1 


as A 


A+ i«l^ ab 


ernents of the Hamiltonian in Eq. 

where £ a b is the elastic strain, \ = y a b is the defor¬ 

mation potential characterizing the linear response of the 
TLS potential to strain, ‘a’ and ‘b’ denote spatial com¬ 
ponents of the strain tensor, and the Einstein summation 
convention for repeated indices is used. Upon applying 
the unitary transformation which diagonalizes the Eq. 
|B1| to the full Hamiltonian including the strain-induced 
perturbation of the asymmetry and averaging over TLS 
orientations, we arrive at Eq. [5] The free Hamiltonian 
for the phonons is given by 


+ 1/2) (B3) 

q,*? 


where the sum counts modes with wavevector and po¬ 
larization q and rj , and bq V and H qi; are the q/y-mode 
annihilation operator and frequency, respectively. The 
decomposition of the strain field into normal modes, here 
in plane waves, reveal the connection with the annihila¬ 
tion operator 5 ql? 


£ab(x) =- ' 52 (q a .fr,b{ q) + <?b^a(q))y ^ 

+ H.c. (B4) 

where V is volume of the system, p is the material density, 
Ujfe(q) is a unit vector for ^-polarized phonons, and H.c. 
stands for Hermitian conjugate. 
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Appendix C: TLS excited state lifetime 


Before we begin our investigation of TLS-induced ef¬ 
fects upon the phonons we outline some of the phonon- 
induced effects upon the TLSs which play a role in our 
analysis. The most important of these effects is that in¬ 
teraction with the phonons leads to a finite TLS upper 
state lifetime. 

For a TLS of asymmetry A and overlap energy Ao 
Fermi’s golden rule gives the excited state decay rate 


1 

r 


Iri -FAq 
2irph A 

<q 'I ' 


coth 


( E 


\2k B T 


(Cl) 


[551 - H^] where the sum over rj counts decay channels cor¬ 
responding with the various phonon polarizations, and T 
is the temperature of the phonon bath. The polarization 
dependence of the deformation potential and the sound 
speed are accounted for in the ry suffices. For fixed energy 
E the lifetime r has the minimum value T m ; n obtained 
from Eq. | C1 1 by taking Ao —► E. 

The lifetime of the jth TLS is derived using Fermi’s 
golden rule. We begin with matrix elements for upward 
transition of the TLS from the ground \(jj) to the excited 
state 1 6j) 


( ej,riqri 1 |Mj : £(rj)a x j \9ji n qri) — 


H 


An 


2 MqnpV 


E 


<-qr/ 


(C2) 


where n ql) is the number of phonon quanta in the q?y- 
mode. The matrix element for downward transitions is 
given by 


since P e + P g = 1 which gives the decay rate for the 
excited state as C g ^ e + C e ^ g written explicitly in Eq. 

EU 


Appendix D: Resonant phonon absorption by TLSs 
1. Derivation of weak field phonon decay rate 


Eq. [6] is derived using the golden rule. The transition 
rates for an increase in the number of quanta in the q? 7 - 
phonon mode from ?i q?) to n qr) + 1 as well as the decay 
rate from n q)) to n qr) — 1 are given by 


C n 


>71 qri + l 




Ej r)P 


(^qr? H - 


C n 


- 1 V p 

7-1 - y zL 


Ojtf nq 2 

9 t? 2 n hflqr,). 


Ej nP 


(Dl) 


The sum over the various TLSs can be performed 
by assuming the validity of the ergodic theorem which 
states that the volume average is equal to the ensem¬ 
ble average of the TLSs in the thermodynamic limit i.e. 

Aoj.-Ej) -a JdAjdA 0 P/A 0 E(A,A 0 ,E) 
where J- represents the summand and recall that P /Ao 
is the TLS density of states. 

Converting the sum over j to an integral over the TLS 
distribution and completing the integrals over A and Ao 
gives time rate of change of n q)/ 


T 1 1 . £ (l"j )^x,j •> n qri ) — 

I h 


nir,\ 


212 q, } pV 


Ao 

~E 


\J n q»; + 1 - 


(C3) 


After averaging over the initial (thermal) phonon state 
and summing over all final states that contribute to the 
two processes the golden rule gives the upper transition 
rate as 


Cg^e — 


7 g EjAl j 


V 


27 Tph 4 


e k B T _ ]_ 


(C4) 


where the Debye phonon density of states has been used 
g v {E) = E 2 /(2 tt 2 H 3 v 3 )V. The rate C e _>. g can simi¬ 
larly be computed and is identical to C g ^ e if the factor 
(exp(i?j/fcsT) —1) _1 is taken to (exp(Ej/k B T) — 1) _1 +1. 

These two rates can be combined to give the time rate 
of change for probability of the TLS to be in its excited 
state 


Pe — Cg—>ePg C e ^ g P e 

= ~ (Cg^ e + C e —>g)P e + C, 


y g->e 


(C5) 


hq T) —A[-P e ( ^q 3 T Pg ^q?/ ] (D2) 

(i) («) (iii) 

where A (P e — P g ) is given by Eq. [6] The under¬ 
braces denote the terms contributing to; (i) stimulated 
phonon emission, (ii) spontaneous phonon emission, and 
(iii) stimulated phonon absorption. From the expression 
above it is clear that when P g —> 1 that the TLSs pre¬ 
dominantly attenuate the acoustic wave. Resonant ab¬ 
sorption is suppressed at either high temperature where 
P e = P g . or high intensities where n q7) 1 and P e —► P g . 
However, the physics of the latter case can only be elu¬ 
cidated by working with the full dynamics governed by 
the Bloch equations. 


2. Strong fields 


At high acoustic intensities perturbation theory is no 
longer adequate to calculate the phonon decay rate. In 
this regime one must solve the coupled Heisenberg equa¬ 
tions of motion for the TLSs and the acoustic field given 
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by 


Cz.i — 


&v.i — 


bq V — 


\ M i ■ £( r iK,j (°3) 

^(•®j + ' £( r j)) a x,j ~ Mj ■ £(r j)cr z ,j (D4) 

~^(^j (b*5) 

—ifiq, ; 6 qr) — — ^ , 9qr],j{^jO'z,j + ^A-OjO'xj) 

(D6) 


where g wJ = jntw e **'**■ 

When the acoustic field is driven strongly at angular 
frequency ui ~ fl qr/ the physics can be greatly simpli¬ 
fied. Under these conditions we focus our attention solely 
on the classical steady-state dynamics of the q^-mode. 
We account for the remaining, thermally populated, 
phonon modes with phenomenological damping terms 
and ‘Langevin’ forces that influence the TLS dynamics. 
Strong driving also allows the use of the rotating wave 
approximation (RWA) that we implement by decompos¬ 
ing the strain field and the x and y Pauli operators into 
positive and negative frequencies i.e. £ = £( + ' + £^- 1 
where oc e =F ' a,t , and similarly a x j = 

where oc e Tlult . The operator a y j is similarly de¬ 

composed and a z j is assumed to be time-independent in 
the steady-state limit. 

After these simplifications the coupled equations of 
motion for the driven phonon amplitude /3q V , and the 
mean values for the TLS operators, (<Jij) = Sij, reduce 
to 


S,,j = 


c(±) 

J y,j 


c(±) 

a x,j 


1 2 

— 7jr{Sz,j — W 0,j) + 


&~\r j)S$+H.c. 


(D7) 


_;L oW 4-— R-9^) — —M y( ± )('r-1 S' ■ 

rp + \ r j)^Z,J 


J_o(±) _ Ip.q(±) 


fiqrj — iQqrjfiqri + Fq V ^ 1 ^0,j 9qr/,j&x.i 


where the TLS phenomenological decay rates l/Tf, quan¬ 
tifying the upper state lifetime, and I/T 2 , characterizing 
the dephasing rate, arise from the interaction with the 
thermal phonon field, and additionally in the latter from 
spectral diffusion of a given TLS’s oscillation frequency 
\{Ej — Dj ■ £(r j)) as the static background strain field 
is modified by spin flips of neighboring TLSs [58]. The 


strain amplitude of the q? 7 -mode £ q ^ is given by the co¬ 


efficient of bq V in Eq. B4 with fe q() replaced by /3 qjJ , M, 




is 7j,(A j/Ej ) 1 and Fq V is the magnitude of an external 
drive at angular frequency oj. The thermal equilibrium 
value of S z j given by Wqj = — tanh(E :/ - /2ksT) plays the 
role of a ‘Langevin’ force in this system of equations by 


ensuring the return to thermal equilibrium in the absence 
of driving. 

In steady-state and noting that the equa¬ 

tion for S z j can be solved in terms of £ q ^ giving 


S,j = 


W 0,j 


1 + 


iM y.j TlT ' 2 |y( + ) 12 

n 2 Isq?) I 




i+sf7f 


(D8) 


where Sj = Ej/h — oj and E j = Ej/h + oj. Eq. |D8| is 
the nonequilibrium steady-state value for the population 
inversion at energy Ej under driving by an acoustic beam 
at angular frequency oj. Eq. |D8| can be used to find the 
steady-state solution for S x j which can then be plugged 
into the equation of motion for the phonon annihilation 
operator to give the TLS-influenced phonon dynamics 


i(^iqrj ca)/f q 


£ 


V,3 


q 2 T x 


2 Hflq^pV 


1 + iT 2 5j 


1 - iT 2 Tij 


S z ,j Pqry 


iAuJan o riJ 


(D9) 


The backreaction of the TLSs on the phonons is charac¬ 
terized by a complex susceptibility — zAtd qr) — |T q ^ that 
modifies the phonon dynamics in two ways: The real part 
results in dissipation of the phonon beam by resonant 
absorption, and the imaginary part induces a frequency 
shift (to be discussed below). 

Assuming the validity of the ergodic hypothesis the 
sum over j can be converted to an integral. By assuming 
that T 2 is insensitive to Ao, and by working in ‘polar’ 
coordinates {E,<j>), i.e A = Ecoscj) and A 0 = E sin (f>, 
the dissipation rate can be expressed as 


pres 

1 qr/ 


P'yffiqy 

hpv 2 


dE- 


T 2 tanh(2? /2ksT) 


l + (f -u>)*T* + 


4 7 


h 2 


— I^| 2 

(DIO) 


where ojT 2 S> 1 has been used, and Tf is taken to be 
approximately given by Eq. |C1| For ujT 2 3> 1, the in¬ 
tegrand is sharply peaked for £ a w. Assuming that 
T 2 , T m i n , and tanh(£?/2/csT) vary little over the range 
E G h{uj — 1 /T 2 ,oj + 1/T 2 ) the integral is given approx¬ 
imately by Eq. [7] where fl ql) = S2, oj = fl, the expres¬ 
sion for the acoustic intensity J = 2pu 3 |£ q ^| 2 has been 
used, and the critical intensity is actually inversely pro¬ 
portional to r m i n as opposed to T). With /(A, A 0 ) and 
Eq. |C1| the density of TLSs as a function of E and r can 
be derived. This distribution is highly peaked near r m ; n 
(see Eq. 3.8 of [42]), and thus we approximate T\ > r m i n 
in Eq. [8] 


Appendix E: Relaxation absorption 

Relaxation absorption results from the modulation of 
the TLS’s energy splitting by a time-dependent strain 
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field. Such an effect was neglected in taking the RWA 
to arrive at Eqs. D7 where we have dropped a S z j- 
dependent drive term in the equation of motion for the 
phonon amplitude. The time dependence of S z j can be 
properly accounted for, and its effect on the phonon dy¬ 
namics can be taken into account, leading to dissipation 
and a frequency shift of the phonons. 


The steady-state amplitude of the time-dependent 
component of S z j that dominantly couples to /3 ql) is given 
by 




dw { 


z ’ j 1-wATi dE, 




(El) 


where SS z j = S z j —wqj. Accounting for this additional 
drive term in the equation of motion for the phonons 
leads to 



u> - 


A 


N , r res 
q?? 




1 \ - 


Di 


Q r 


dw, 


o,j 


V ^ 4pr>2 1 - i w Ti dEj 




-?Aa) rel -ir rel 

2 i qrj 


(E2) 


Just as in the case of resonant absorption, the effect of re¬ 
laxation absorption is quantified by the complex suscep¬ 
tibility -iAw” 1 — 2 TqjJ. Here we focus on the relaxation- 
induced damping coefficient 


prel 

1 q^? 


2 D^ j dwgj 

V ^ 4pz| 1 + oAT 2 dEj ■ 


(E3) 


To evaluate T^j we employ the ergodic hypothesis and 
note an important observation about the relative magni¬ 
tudes of the probing frequency and the inversion decay 
in our experiment. Assuming that T\ is given by the up¬ 
per state lifetime of the TLSs then it takes a minimum 
value calculated by taking Ao —> E in Eq. |C1| Fur¬ 
thermore, the factor dw 0 /dE exponentially suppresses 
contributions to the integrals from E bigger than ksT. 
Hence, for the frequencies of interest in our experiment 
it can be shown that wT) 3> 1 over the integration range 
contributing to T^. Thus, a Taylor expansion of the 
integrand of T^ is justified for small 1/wXj and results 
in 


/ y ^ 

jnrel _ ly 

qv pv*k B T 


dA 


dAo 


A 2 3 4 5 6 7 8 9 10 11 1 1 

A 0 E 2 t cosh 2 


(E4) 


where we’ve taken f l qv /ui = 1 and where we’ve replaced 
Tj with r. One arrives at Eq. 11 after performing the 
integrals. 
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